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' Electronic structure calculations were carried out on the compound GdBaCo2 05. 5. The electronic 

, structure variation with a change in the spin state of the Co''"'' ion in an octahedral environment 

■ has been studied. All the insulating phases are described and possible scenarios for the metallic 

^Nj ' ones are presented. Orbital ordering is shown to take place and the electronic structure leading 

» I to it is determined. The Ising-like anisotropy shown experimentally can be predicted. Also, big 

. unquenched orbital angular momenta are calculated and their origin is described. 
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PACS numbers: 71.20.-b, 71.30.+h, 75.50.Pp 



I. INTRODUCTION 



Co oxides have been the objective of deep research 
during the last few years. As any correlated oxide, 
they present an interplay between the electronic struc- 
ture, magnetic properties and geometric structure that 
makes them very interesting from both an experimental 
and a theoretical point of view. They have received a 
strong attention very recently, mainly after the discov- 
ery of both superconductivity^ and magnetoresistance^ 
in them because of the vast technological applications of 
those phenomena, but also because of the magnetic tran- 
sitions accompanied of structural transitions, that are 
not yet fully understood. GdBaCo205.5 presents mag- 
netic transitions at about 75 and 220 K and a metal- 
insulator transition at some 360 K?^ The magnetization 
measurements at the different phases have not been con- 
clusively ascribed to a particular spin state and orbital 
configuration. There are no measurements available of 
the orbital magnetic moments of the Co ions in the 
different phases and spin states. This rich variety of 
phases including spin state transitions occurs in Co^+ 
compounds, where the low-spin state is not stable at all 
temperatures because of the small magnitude of the crys- 
tal field, comparable to the intra-atomic exchange energy. 
This is the case of GdBaCo205.5, that also presents gi- 
ant magnetoresistance^ associated with a metal-insulator 
transition. Very recently, a spin blockade phenomenon 
has been found to explain its conduction properties'*^ 
Some work has also been carried out from a theoreti- 
cal standpoint iSiL2iS The goal of this paper is to make 
use of electronic structure calculations in order to ex- 
plain the magnetic and electronic structure of the com- 
pound and their variations with respect to the magnetic 
and spin state transitions that occur in it. For doing so, 
in section ^] we will describe the structure of the mate- 
rial and the method of calculation and in section IIIII we 
will present our calculations of the magnetic and elec- 
tronic structure, dealing with the spin state, magnetic 
and metal-insulator transitions, and also with the orbital 
ordering phenomenon relating all of these to the macro- 
scopic observations establishing a consistent picture for 
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FIG. 1: (Color online) Structure of GdBaCo205.5. Grey 
atoms (Col and Co2) are Copj,,-, dark atoms (Co3 and Co4) 
are Cooct and the small light ones are the oxygen atoms. Ob- 
serve the oxygen vacants along the c-axis leading to a square 
pyramidal environment. 



the observed properties of the material. 



II. STRUCTURE AND COMPUTATIONAL 
DETAILS 

GdBaCo205.5 presents two crystallographically differ- 
ent surroundings for the two inequivalent Co ions in the 
unit cell. One of them is in a square pyramidal (Copj,^) 
and the other in an octahedral environment (Cooct). 
There is one oxygen vacant per unit cell, that produces 
the square pyramidal environment for the Co in the crys- 
tallographic position 2r. The octahedron surrounding 
Cooct is elongated along the a-axis (see Fig. ^for the 
naming of the crystallographic directions). The unit cell 
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is orthorhombic with space group Pmmm. The value of 
the lattice parameters and atomic positions were taken 
from Respaud et al}° {a= 3.87738 A, b= 7.53487 A and 
c= 7.8269 A), but the oxygen positions were recalculated 
performing a full structural relaxation. It is well estab- 
lished that the oxygen vacants are in the crystallographic 
position Ig (1/2, 1/2, l/2)m^ For our calculations we 
broke the symmetry along the b-axis, so that 4 inequiv- 
alent Co atoms enter the cell (let us call Col and Co2 
the two inequivalent atoms in the pyramidal position and 
Co3 and Co4 the two inequivalent atoms in the octahe- 
dral environment, as can be seen in Fig. U}. Further 
symmetry breakings were produced along the c- and a- 
axis to have 8 Co atoms in the supercell: the former 
case with the intention of studying possible spin state 
orderings in the yz plane and the latter for analyzing the 
magnetic couplings between Copyr layers mediated by a 
plane of Cooct atoms. 

TABLE I: Crystallographic positions at T= that result from 
our structure optimization. 



Atom 


Crystallographic 
position 


Coordinates 


Ba 


2cr 


(0,0,0.2503) 


Gd 


2p 


(0,0.5,0.2292) 


COpyT- 


2r 


(0.5,0.2563,0.5) 


COoct 


2q 


(0.5,0.2487,0) 


01 


la 


(0.5,0,0) 


02 


le 


(0.5,0,0.5) 


03 


Ic 


(0.5,0.5,0) 


04 


2s 


(0,0.2731,0) 


05 


2t 


(0,0.3102,0.5) 


06 


4u 


(0.5,0.2947,0.2633) 



We present here full-potential, all-electron, electronic 
structure calculations based on the density functional 
theory (DFT) utihzing the APW-|-lo method" per- 
formed with the WIEN2k software liSii^ For the structure 
minimization, the results of which are summarized in Ta- 
ble ^ we used the GGA (generalized gradient approxima- 
tion) in the Perdew-Burke-Ernzerhoff (PBE) scheme.iA 
The geometry optimization was carried out minimizing 
the forces in the atoms and the total energy of the system. 
These results are independent of the magnetic ordering, 
forces are minimized both in the case of ferromagnetic 
(FM) or antiferromagnetic (AF) ordering. Also, a force 
minimization procedure was carried out in a supercell 
containing 8 inequivalent Co atoms. Minimal changes of 
the atomic positions compared to those reported in Table 
Uturn out as a result. For the electronic structure calcula- 
tions, we included the strong correlation effects by means 
of the LDA-I-U schemc^'"^ in the so-called "fully-localized 
limit" The values of U chosen were within the range of 
3 to 6 eV for all the spin states (U depends on the spin 
state of the atom) but all the results presented here are 
fully consistent for U within that range (the quantitative 
results of energy differences depend on the choice of U by 



some 25% for values of U within that range, the values 
presented will be given for U= 4.5 eV). Spin-orbit effects 
have been introduced in a second variational way using 
the scalar relativistic approximation.^^ The parameters 
of our calculations depend on the type of calculation but 
for any of them we converged with respect to the k-mesh 
and to Kmt^maxi up to 1000 k-points (256 in the irre- 
ducible Brillouin zone) and up to RmtKmaa;= 7. Local 
orbitals were added for a bigger flexibility in dealing with 
the semicore states. Mufhn-tin radii chosen were the fol- 
lowing: 2.3 a.u. for both Ba and Gd, 1.90 a.u. for Co 
and 1.59 a.u. for O. 



III. ELECTRONIC AND MAGNETIC 
STRUCTURE 

From a purely ionic point of view, we expect both Co^"*" 
ions in a d^ configuration, that in an octahedral environ- 
ment could lead to either a low spin (LS) state (t2ge"; 
nonmagnetic S=0), an intermediate spin (IS) state (t^gci; 
S=l) or even a high spin (HS) state (tl^e^; S=2). Spin 
state transitions play an important role in the temper- 
ature evolution of the electronic, structural and mag- 
netic properties of transition metal oxides, specially in 
the case of Co'^"'" compounds, such as the famous example 
of LaCoOa.^^ In the case of GdBaCo205.5, it is not yet 
clearly understood the relationship between a spin state 
transition and the transition from AF to FM behavioria 
(neither with the metal- insulator transition). We have 
studied in this paper how the electronic structure evolves 
with spin states and our calculations always relax to 
a scenario with Copj,^ in the IS state, not being pos- 
sible the stabilization of a different spin state for that 
atom. Usually, the pryramidal sites of GdBaCo205.5 
have been ascribed to the IS statej^LSSi even crystal field 
arguments confirm it as being very stabled Very recently, 
a spin blockade phenomenon was found to explain the 
conduction properties of the material,"* being its origin 
explainable assuming an IS state in Copyr atoms^ How- 
ever, some works evidence that a HS state is likely to 
occur for a Co'^"'' ion in a square pyramidal environ- 
ment in some cases, as has been evidenced by neutron 
diffractioufS^iSi spectroscopy measurements jiS^ and also 
theoreticall}i^i2S42i for different materials containing sim- 
ilar environments for a Co''+ ion. Even for GdBaCo205.5, 
calculationsSii and theoretical results^ show a HS state 
as the most stable for Copyr, contrary to our results. 
However, all the solutions we tried converge to an IS 
state of Gopyr indicating an enormous stability of such a 
state. The different results obtained with similar meth- 
ods could be due to our structure optimization, that leads 
to changes of up to 2% in interatomic distances and bond 
angles with respect to the experimental room tempera- 
ture data. It is well-known that variations in the geome- 
try of the cation environment could lead to the stabiliza- 
tion of a different spin stateS^ Also, the LDA-I-U method 
is strongly dependent on the starting point and several 
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solutions can be obtained depending on the starting den- 
sity chosen. On the other hand, for Cooct; different spin 
states could be converged: a LS state, an IS state and a 
spin state order with half the Gooct atoms in a HS state 
and the other half in a LS state. A HS state of Cooct 
could not be converged with the structural data we uti- 
lized. It has been shown experimentally that a structural 
transitioniS occurs at the metal-insulator transition and 
this has been related to the appearance of a HS state. 
Such a HS state would stabilize through a strong spin- 
lattice coupling that our zero temperature DFT calcula- 
tions cannot describe. 



A. Magnetic ordering and magnetic transitions 

We have studied several magnetic and spin configu- 
rations within the LDA-I-U approximation. Both FM 
and AF couplings along the three crystallographic di- 
rections have been considered. Our calculations yield a 
FM coupling along the b- and c-axis as the most stable 
one for any spin state of the COoct layer. Within each of 
the layers, an AF ordering cannot be stabilized. It has 
been argued^ that orbital ordering produces the stabi- 
lization of a FM ordering within the plane, and that is 
what we observe (see Section Fill B|) . For the case of a 
LS state of the Cooct atoms, a FM and an AF ordering 
of Copyr layers mediated by a nonmagnetic Cooct layer 
were calculated. AF ordering is more stable by some 10 
meV/Co (U-dependent value but consistent for a wide 
set of U choices), not a strong magnetic coupling. Both 
IS-LS solutions (IS stands for the spin state of Copyr and 
LS for the spin state of the Cooct atoms), consisting of 
planes coupled AF or FM, lead to an insulating behavior 
(see Fig. We observe that the IS-LS configuration 
leads to a zero- gap insulator, in accordance with experi- 
ments observing a narrow gapi^S This is the only solution 
that could explain the properties of the AFI (antiferro- 
magnetic insulator) phase found below 70 K. The weak 
magnetic coupling between the layers breaks up at that 
temperature entering the paramagnetic phase above it. 

From our calculations, we interpret the AF-FM tran- 
sition that occurs at about 220 K as a change in the 
magnetic coupling of FM Copj,^ layers mediated by non- 
magnetic Cooct planes, because such a FM ordering is 
the only one leading to a FMI (ferromagnetic insulator) 
phase. This is consistent with the experimental observa- 
tions in Ref. 

The magnetic coupling between Cooct and Copj,^ lay- 
ers when the former acquires an IS state is FM. The en- 
ergy difference between a FM and an AF configuration 
is about 40 meV/Co (U-dependent value), a strong mag- 
netic coupling. We could also converge a FM solution 
where half of the Cooct atoms are in a HS state and the 
other half in a LS state, a spin state ordering situation, 
similar to what has been observed for NdBaCo205.g^ 
(in this case it has been predicted to be between IS 
and LS states) or predicted theoretically for the whole 
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FIG. 2; Spin-up and spin-down total DOS plots for the ma- 
terial in the insulating phases, when the Cooct atoms are in 
a LS state. The upper panel shows an AF coupling between 
Copyr layers and the lower panel shows a FM situation. 



serieSiS^ (but in that case the authors assume that Gopyr 
posses a HS). These two solutions we have computed lead 
to a half- metallic state (see Fig. ISJ. In principle, any 
of these could be the magnetic ordering in the metallic 
phase above the metal- insulator transition. But, since a 
structural deformation has been observed experimentally 
at the metal-insulator transitioniS and this has been as- 
cribed to the onset of a HS state, we should not rule 
out the possibility of a HS Cooct layer above the metal- 
insulator transition. The structural data we utilized (cal- 
culated at T= 0) does not allow to stabilize that config- 
uration. The calculation of such a phase by DFT tech- 
niques would require a precise knowledge of the geometry 
above 360 K. The recent magnetization measurements^ 
observing a total magnetic moment of 1 (Ib/ Co above the 
metal-insulator transition would rule out this possibility. 
As we will see below, orbital angular momenta need to be 
taken into account for explaining the total magnetization 
measurements . 



B. Orbital ordering 

Co^+ ions in an IS state have an orbital degree of free- 
dom, that could lead to the existence of orbital order- 
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FIG. 3; Spin-up and spin-down total DOS plots for the ma- 
terial in the possible metallic phases considered in the text. 
The upper panel shows a spin state ordering HS-LS in the 
Cooct layers and the lower panel shows an IS state for all the 
Cooct atoms. 



ing along the b and c-axesr^S^Si The octahedral environ- 
ment is elongated along the crystallographic c-axis (we 
will consider it here as our z-axis for the representation 
of the d-levels). For the pyramidal environment, the z- 
axis will be considered along the direction of the oxygen 
vacant (crystallographic b-axis). For Copyr, due to the 
small crystal field splitting a Co'^'^ ion suffers, an orbital 
degeneracy is likely to occur as well (depending on the 
intra-atomic exchange energy). 

Let us first begin with the IS-LS situation, the ground 
state. Figs. 0] and |S1 show the density of states (DOS) 
of the d-levels for the IS state of the Copyr atoms. Their 
electronic structure does not vary much when the spin 
state of Cooct changes. According to our calculations, 



the IS state of Copyr is fairly stable and remains un- 
changed when increasing temperature, so our description 
of this state is valid for any spin state of the Cooct atoms. 
We can observe in the curves a different electronic struc- 
ture in the crystallographically equivalent Co atoms, in- 
dicating the existence of orbital ordering along the 6-axis 
(the crystallographic axis that joins the two inequivalent 
Copyr (and also Cooct) atoms we have considered in our 
calculations (see Fig. 0. In Ref. ^3, this orbital state 
is predicted to be dxz i i dyz for DyBaCo205.5, but we 
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FIG. 4: Spin-up and spin-down partial density of states of 
Col (in a square pyramidal environment). It can be observed 
that it is in an IS state and the hole in the t2g multiplet is in 
the dyz orbital. 



obtain as the ground state an orbital ordering between 
real orbitals (d^z and dyz), being this configuration some 
20 meV/Co more stable than the state with equal orbital 
structure along the 6-axis (even though the magnitude of 
the energy difference is dependent on the value of U cho- 
sen, the important point is that it certainly is more stable 
independently of U), a value comparable to the energy 
involved in magnetic ordering. Similar calculations indi- 
cate that the same orbital ordering is stable as well along 
the c-axis. Couplings within the bc-plane are equivalent 
both orbitally and magnetically. However, the location 
of the electron that is promoted into the eg doublet is 
more difficult to ascribe, since it is not in a pure d322_r2 
nor in the dj.2_y2 orbital, as can be observed in Figs. ^ 
and [51 We also decomposed the DOS curves in a set of 
orbitals using the complex functions d^z^-r^ i i d^^-y-^ 
(consistent with having the same environment for both 
Co atoms, since these orbitals have cubic symmetry) but 
the electron does not fit into those orbitals alternately 
neither. Some more complicated orbital structure occurs 
that we could not discern with our calculations. We can 
only ascertain that the resultant states must have some 
component of each orbital, probably a complex combi- 
nation since an unquenched orbital angular momentum 
occurs as we will see below. Also, the Bg bands (this is 
valid for Cooct in an IS state as well) contain more than 
one electron, that probably comes from a very strong 
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FIG. 5: Spin-up and spin-down partial density of states of 
Co2 (in a square pyramidal environment). It can be observed 
that it is in an IS state and the hole in the t2g multiplet is in 
the dxz orbital. 



FIG. 6: Spin-up and spin-down partial density of states of Co3 
(in an octahedral environment) for the IS-IS configuration. It 
can be observed that it is in an IS state and the hole in the 
t2g multiplet is in the d^cz orbital. 



hybridization of the band with the surrounding O p 
bands. This could imply the formation of a Co^+L state 
(L being a hgand hole). This strong hybridization had 
been already found in previous worksfSiS but for a HS 
state of Copyr- 

For the case of Cooct, we present in Figs. El and [3 the 
d-DOS for the IS state. The plots of the LS state are not 
introduced because the electronic structure is basically a 
^2g configuration in a nearly octahedral environment, 
the t2g levels fully occupied and the eg levels fully unoc- 
cupied. In the IS state, one of the electrons in the lower 
lying multiplet is promoted to the e^ levels. There oc- 
curs an orbital ordering phenomenon as well for Cooct in 
the IS state along the 6-axis. The DOS of two contiguous 
Cooct atoms along that direction are shown in Figs. |^and 
|7| where it can be seen that the hole left in the t2g mul- 
tiplet is alternatively at the d^y and d^.^ orbitals. This 
configuration is some 10 meV/Co more stable than the 
non-orbitally ordered scenario (that energy value might 
vary slightly with the value of U chosen, but it is con- 
sistent for a wide range of U) . The electron that is pro- 
moted in the Cg multiplet places in alternate orbitals, but 
from our calculations it is difficult to describe the com- 
bination of c?322_r2 and dx2_y2 orbitals that forms the 
doublet. We know it is not the pure real orbitals as we 
can see in Figs. |51 and [7| but they are not the complex 
orbitals dy,^-^^^'^ ± i d,j.2_y2 neither. We tried that possi- 



bility decomposing our DOS in a basis containing those 
two orbitals but such an ordering was not found. Some 
other admixture of these orbitals is happening due to the 
local environment of the atom but from our calculations 
this could not be sorted out. We encounter here a similar 
situation to the case of Copyr, the fact that there exist 
some unquenched orbital angular momenta rules out the 
possibility of having real orbitals, the electron must be 
located in a state formed by some complex combination 
of orbitals. The e^ band is also in this case strongly hy- 
bridized, possibly forming a Co^+L state as for Copyr- In 
an IS-IS situation, the orbital ordering was found to be 
stable both along the b- and c-axis, consistent with the 
picture established in Ref. |^ 

In a spin state ordering situation, the Cooct atom in 
a LS state acquires some magnetic moment due to the 
surrounding HS ions and the Cg bands become occupied 
forming a very wide band that crosses the Fermi level 
leading to the half-metallic behavior observed in Fig. |3| 



C. Spin-orbit study: unquenched orbital angular 
momenta and magnetic anisotropy 

We have also carried out calculations including spin- 
orbit coupling to elucidate how the magnetic anisotropy 
changes for the different spin states and magnetic cou- 
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FIG. 7: Spin-up and spin-down partial density of states of Co4 
(in an octahedral environment) for the IS-IS configuration. It 
can be observed that it is in an IS state and the hole in the 
t2g multiplet is in the d^y orbital. 



plings and to evaluate the orbital angular momenta 
(0AM) of the different Co atoms in the structure, that 
had been predicted to have some unquenched OAM^i but 
whose values have never been measured. Only values of 
total magnetization are available in the literature, ^^'^i^ 
showing some discrepancies between them. 

TABLE II: Values of the orbital angular momenta of the four 
Co atoms in the structure in the three spin state configura- 
tions studied. The values are given in fj,B and are considered 
inside the muffin-tin sphere. These values are calculated for 
the magnetization along the experimental easy axis. 



Atom 


IS-IS (/is) 


IS-LS (hb) 


IS-(HS-LS)(/iB) 


Col (pyr) 


0.30 


0.28 


0.27 


Co2 (pyr) 


0.02 


0.01 


0.01 


Co3 (oct) 


0.01 


-0.05 


0.03 (LS) 


Co4 (oct) 


0.24 


-0.07 


0.19 (HS) 



A summary of the 0AM calculated is presented in Ta- 
ble m These data are calculated assuming the magne- 
tization goes along the experimental easy axis (c-axis). 
The values are given inside the mufhn-tin spheres, so the 
real values of the moments are, hence, somewhat bigger 
(about 10-30%, but this is difhcult to estimate). They 
also are U-dependent by up to 25% for other combina- 
tions of U parameters within a reasonable range. Large 



values up to 0.6 /is for some of the Co atoms would be 
consistent with our calculations. Experimental measure- 
ments are needed to confirm this observation. We per- 
formed these calculations for the three different spin state 
configurations we could converge for the Cooct atoms 
{Copyr is always in an IS state): a LS, an IS state and 
a spin state ordering with half the Cooct atoms in a HS 
state and the other half in a LS state. 

Irrespective of the Cooct spin state, we observe that 
Copyr atoms have very similar angular momenta (one of 
them being about 0.3 fiB inside the muffin-tin sphere). 
The electronic structure of these atoms barely changes 
when a spin state transition occurs to Cooc*. This atom 
develops a big angular momentum in one of the sites 
when it is in an IS state. It is worth noting that the big 
orbital angular momentum occurs when the hole in the 
t2g multiplet is left in the dxz orbital, both in the case of 
Copyr and Cooct in an IS state, whereas the configuration 
with a hole in a different t2g orbital produces a negligible 
moment. Small values are predicted also for Cooct in a 
LS sate. These are oppositely oriented to the magnetic 
moment of the Copj,^ ions when all the Cooct atoms are in 
a LS state and parallel to them when the layer includes 
Cooct atoms in a HS state. These also have a large value 
of the 0AM, but smaller than those of the IS state. 

We also performed calculations of the magnetic 
anisotropy energies and our results agree with the ex- 
perimental data availableiSS We obtain a strong Ising- 
like behavior for the AFI phase (IS-LS with IS Copyr 
planes coupled AF via nonmagnetic Cooct layers). In 
that case, the c-axis is the easy axis by some 150 K/Co, 
a big anisotropy energy, in agreement with the experi- 
mental value (some 80-100 K/Co).^^ For the IS-LS FMI 
phase, we find an easy axis along the c-axis, being the 
b-direction almost degenerate (3 K/Co harder) and the 
a-direction harder by some 40 K/Co. No trace of Ising- 
like behavior is found for this phase. This disagrees with 
the experimental findings because our zero-temperature 
DFT calculations cannot describe appropriately the mag- 
netic anisotropy properties at 260 K, where they have 
been measured. A similar situation occurs for the IS- 
IS FM state, the b-axis is the easy axis being the a and 
c directions harder by some 30 K/Co. In these cases it 
could be possible that the domain formation found in 
Ref. could have an influence on the value of the mag- 
netic anisotropy constants that is not considered in the 
calculation of the magnetocrystalline anisotropy energy. 



IV. CONCLUSIONS 

In this paper we have presented a study of the elec- 
tronic structure and magnetic properties of the Co oxide 
GdBaCo205.5 using ab initio calculations considering an 
all-electron, full-potential APW-l-lo method. The points 
we tried to address were the spin state transitions, the 
orbital ordering and spin-orbit effects (unquenched or- 
bital angular momenta and magnetic anisotropy proper- 
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ties). From our calculations we confirm the stability of 
the IS state for Copyr or the conduction properties of 
the IS-LS state as a narrow gap insulator. The transi- 
tion to a metallic phase must be accompanied of a spin 
state transition of Coqc*. It is necessary to know the 
precise geometry in the metallic phase to determine its 
spin and magnetic configuration. We have observed the 
material is in an orbitally-ordered state both at low tem- 
peratures in the IS-LS state and also in a hypothetical 
IS-IS state, and we described the orbitals involved in the 
phenomenon. We also predict big unquenched orbital 
angular momenta for some of the Co atoms in the struc- 
ture, a fact that had been predicted but whose values 
have never been measured. We estimate their values and 
explain their origin. Also, the magnetic anisotropy was 



studied and agrees with the experimental observation of 
a strong Ising-like behavior. In summary, we establish a 
consistent picture for the properties of the material by 
confirming some experimental evidences and predicting 
new results based on first principles calculations. 
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